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In this talk we discuss the microscopic limit of QCD at nonzero chemical po- 
tential. In this domain, where the QCD partition function is under complete 
analytical control, we uncover an entirely new link between the spectral density 
of the Dirac operator and the chiral condensate: violent complex oscillations 
on the microscopic scale give rise to the discontinuity of the chiral condensate 
at zero quark mass. We first establish this relation exactly within a random 
matrix framework and then analyze the importance of the individual modes 
by Fourier analysis. 
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1. Introduction 

Many firm results about QCD are based on the fact that the Euclidean 
partition function has a real and positive measure. Prime examples are 
QCD inequalities and the evaluation of the lattice QCD partition function 
by Monte Carlo simulations. For this reason much of the intuition gained 
about nonperturbative QCD is based on a probabilistic interpretation. If 
we introduce a chemical potential in order to favor the presence of quarks 
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over anti-quarks the measure in the Euclidean partition function is no longer 
positive definite. In this lecture we will analyze the link between the spectral 
density of the QCD Dirac operator and the order parameter for spontaneous 
chiral symmetry breaking. We will see that the intuition based on zero 
chemical potential results fails completely at nonzero chemical potential. 

At zero chemical potential the eigenvalues of the anti-Hermitean Dirac 
operator are located on the imaginary axis. The chiral condensate can be 
thought of as the electric field at the position on the real axis given by the 
quark mass, to, created by positive charges located at the purely imaginary 
eigenvalues [1]. The chiral condensate is the jump in the electric field as 
the quark mass passes through zero. Clearly such a jump will occur if the 
density of the eigenvalues/charges near the origin of the imaginary axis 
scales with the space-time volume [2]. 

At nonzero chemical potential the anti-Hermiticity of the Dirac opera- 
tor is lost and the eigenvalues are located homogeneously in a strip parallel 
to the imaginary axis with width proportional to fj 2 (for sufficiently large 
^i). The constant of proportionality is such that, at sufficiently low tem- 
peratures, the quark mass hits the eigenvalue strip for /i = [3-5]. 
From our intuition at zero chemical potential we are lead to the conclusion 
that at any (however small) nonzero value of the chemical potential the 
chiral condensate has no discontinuity when the quark mass crosses the 
imaginary axis. This is in sharp contrast to the expected phase diagram, 
see Fig. 1. (Notice however, that in phase quenched QCD the chiral con- 
densate rotates into a pion condensate so that chiral symmetry remains 
broken spontaneously in spite of a vanishing chiral condensate.) Increasing 
the temperature tends to decrease the width of the strip of eigenvalues so 
that the value of the chemical potential for which the quark mass is inside 
the domain of the eigenvalues increases as well (see Fig. 1 for illustration). 

At nonzero chemical potential the eigenvalue density of the Dirac oper- 
ator is defined as 



Since the measure includes the fermion determinants there is no guarantee 
that the density is real or positive. To completely destroy the electrostatic 
picture it turns out that the eigenvalue density depends strongly on the 
quark mass [7,8]. Hence we can not interpret the quark mass as a test charge 
which, in a harmless manner, measures the chiral condensate/electric field. 
To study the eigenvalue density we turn to the microscopic scale of QCD [9] 
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Fig. 1. Schematic phase diagram of QCD and region where the quark mass is inside the 
support of the eigenvalue density. The dashed curve indicates the chiral transition be it 
first order or crossover. To the right of the full line the quark mass is inside the support 
of the Dirac spectrum. The solid curve does not indicate a phase transition in QCD. 



where the dimensionless combinations 



are kept fixed as the volume V is taken to infinity (S is the chiral condensate 
and F-ji the pion decay constant). In this limit the Compton wavelength of 
the pions is much larger than the linear extent of the volume and the zero 
momentum mode of the pions dominates the low energy effective theory 
of QCD [10,11] resulting in a partition function that is given by a group 
integral. Here we will access the microscopic limit of QCD by analyzing 
the chiral random matrix model introduced previously [7]. We will show 
that the unqucnched eigenvalue density is complex and oscillating when 
/i > m n /2 and that these oscillations are responsible for the discontinuity 
of the chiral condensate [12]. The complex eigenvalue density of the Dirac 
operator for QCD at nonzero chemical potential was also studied for gauge 
fields given by a liquid of instantons [6] and strong fluctuations in the 
spectral density could be identified. 



771 = raYjV and fi 2 = /i 2 F 2 V 



(2) 
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2. The random matrix model 

The random matrix partition function with Nf quark flavors with mass m 
and chemical potential /i is defined as [7] 

Z^ s {m;ii)= j d<Pd^> w G ($)w G (^)det N f {V(n) +m), (3) 

where the non-Hermitian Dirac operator is given by 

^ = (»t + M *t o )■ (4) 

The (N + v) x N matrices $ and ^ are complex, v is the topological index, 
and wq is a Gaussian weight function 

wgPO = exp(-NTiX^X) . (5) 

In this random matrix model, the microscopic limit is defined as the limit 
N — > oo where 

m = 2iVm and /i 2 = 2N^i 2 (6) 

are kept fixed. The random matrix partition function becomes identical to 
the microscopic QCD partition function provided that we identify, see for 
example [8], 

m = 2Nm -> m^V, (7) 
A 2 = 27V/i 2 - M 2 F 2 V. 

The main advantage of this matrix model as compared to the one originally 
introduced in [13] is that an eigenvalue representation is known [7] 

r N 

^(m;/i)~m v ^ / J] d 2 z k V N '({ Zi }, K},m; M ), (8) 
^ fc=1 



where the joint probability distribution reads 

1 N 

r N '({ Zl },{z:}, m ^) ^ — \A N ({zf})\ 2 Hw(z k ,z* k ;n)(m 2 -zi) N f. 

V fc=i 
The Vandermonde determinant is defined as 

iV 

A^ 2 })^ LI ( 9 ) 
i>j=l 

and the weight function includes a modified Bessel function, 

W (z k , 4 ;») = \z k \^K v (^^W?) e-^^'K (10) 
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Given the eigenvalue representation we can employ the powerful method of 
orthogonal polynomials in the complex plane developed in [14-19]. 

2.1. Orthogonal polynomials, partition function and chiral 
condensate 

The complex orthogonal polynomials corresponding to the weight function 
(10) can be expressed in terms of complex Laguerre polynomials [7] 

»(,;,) = (^)V <") 

They satisfy the orthogonality relations 

d 2 z w(z,z*;n) p k (z;n) Pi(z; u)* = S k i r£ , (12) 

with the norm given by 

„ _ (l+^ 2 ) 2fc +- fc! {k + v)\ 

The partition function for one fermion can be expressed in terms orthogonal 
polynomials as 

Z^ /_1 (m;^) = m v p N (m;u). (14) 
Using this partition function (14) we immediately find the chiral condensate 
^N f =i, \ d 7 N f =i, > dp N {m)/dm v 

^ N ( m ) = -7- lo g^Ar ( m ;M = r~\ — + — ■ ( 15 ) 

dm PN{m) m 

Since we have not yet taken the microscopic limit the chiral condensate 
depends on p. 



3. The chiral condensate from the eigenvalue density using 
orthogonal polynomials 

The quenched eigenvalue density given by 

n s =o, * s / * , Pk{z*)pk{z) 

k=o Tk 

is real and positive. The unquenched spectral density [7] which can be 
expressed as 

P T\z, z\m- a) = W (z, z*\ a) £ Mz')Wz)-PN{z)p k {m)lp»{m)) 



(17) 
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is not invariant under complex conjugation. The asymmetric second term 
causes complex valued oscillations, see figure 2. 

Here we will compute the chiral condensate (15) starting from the com- 
plex and oscillating eigenvalue density given in (17), 



E^' 1 (m)= f d 2 z 2m p N N ! \z,z*,m;n). 
Jc z " l 



(18) 



We will show that all integrals can be carried out using the orthogonal- 
ity condition (12). The first step is to insert the above expression for the 
unquenched density 



jf f=1 (m) = [ d 2 z 2mW ^ Pk( z *)(Pk{z) ~PN(z)p k (m)/p N (m)) 
N Jc 7 - 2 - m2 f^y r k 

We now use the fact that 

f d2z w(z 7 z*;fi)p k (z*)p k (z) = r d2 ^ w(z,z*;fi)p k (z*)p k (m) 
Jc 7 - 2 - m2 Jc z 2 -m 2 ' 

which is easily proved by first rewriting (the c k s are independent of z) 

p k (z) = (z 2 - m 2 )[p fe _i(z) + c k - 2 Pk-2(z) + ■■■} +Pk(m) (21) 
and then using orthogonality, and find 

^JV/=i/ \ f ,2 2mw Pk{z*)p k {m) 1 

n M = / d z— j 2^ T^iPN {m )~PN (z (22 

Jc z m t~f Tk PN(m) 

For PN{m) — pn{z) we now insert (this is just Eq. (21) rearranged a bit) 

p N (m) -pn(z) = -{z 2 - m 2 )[p N -i(z) + Cj- 2 p N -2(z) + ...] (23) 

and get 

r N f =1 / \ f a2. o Pk{z*)pk{m) 

T, N ' (m) = - / d z z2mw) ——[p N _ 1 (z)+c N ^ 2 PN-2{z) + ...}. 

Jc £r* r kPN {m) 

(24) 

Using orthogonality we find (note that the ^-dependent norms drop out) 

S^ /_1 (to) = 2m [pjv-i(m) + c N - 2 PN-2(m) + ...}. (25) 
PN{m) 
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Re[z] 

Fig. 2. The eigenvalue density of the random matrix model (3) for TV = 20, fi = 0.8 and 
v = 0. Top: The quark mass is well outside the support of the spectrum at m = 2.0, and 
the eigenvalue density is real and positive. Bottom: The quark mass is now inside the 
eigenvalue cloud at m = 0.6, and oscillations starting at z = ±m extend outward. Only 
the real part of the spectral density is shown - the imaginary part is nonzero inside the 
oscillating region. 
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Now, 



dp N (m) _ p N (z)-p N (m) 



dm z^m z — m 

,. o Pn(z) -J)jv(m) 
= lim 2m ^ 5 , 

z^m z — m 

= lim 2m[p N -i(z) + c N - 2 PN-2{z) + ■■■}, 

z—>m 

= 2m[p N ^ 1 (m) + c N ^ 2 PN-2(m) + ...}, (26) 



so that 



S rv)=^4 (27) 



This is the desired relation (15) up to the term v/m which is due to the 
contribution of v exact zero modes that were not included in the spectral 
density. The result of this section originally appeared in [12]. 

Note that the proof relies entirely on orthogonality properties of poly- 
nomials. The detailed cancellations leading to the simple form of the chiral 
condensate are in this way linked to the orthogonality properties in the 
complex plane. 

The above proof was given for finite N and is of course also valid in 
the universal microscopic limit. As shown in [12], the proof can also be 
carried through in the microscopic limit without relying on the finite TV 
chiral random matrix model. That is, one can start from the microscopic 
limit of the spectral density (17), which is [7,8] 

p^ =1 (z,r>;A) = ^K v (M!) (28) 



(***) . 



x ( [ dtt e- 2 ^ 2 * 2 ' Uzt)IJz*t) - ^fer t dt t e" 2 ^*' ' IJmt)I v 
\Jo tv{m) Jo 

and compute the chiral condensate using complex contour integration tech- 
niques. In the microscopic limit the cancellations are even more dramatic 
than at finite N. The entire p dependence of the eigenvalue density vanishes 
upon integration and leaves us with a fi independent chiral condensate 

s <*> = tf| «> 

Rather than repeating the proof of [12] we now focus on the Fourier com- 
ponents of the eigenvalue density that are responsible for chiral symmetry 
breaking. 
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Note that the first term in (28) is the quenched eigenvalue density [20] 
and the effect of unquenching is in the second term. This motivates the 
notation 

P%U=P ( q ) -P ( u\ (30) 
which we shall make use of below. 

4. The Fourier transform of the eigenvalue density 

To further expose the nature of the cancellations leading to a p independent 
chiral condensate we now consider the Fourier transform of the eigenvalue 
density. The motivation to do so is that the discontinuity in the chiral 
condensate is due to the strongly oscillating part of the spectral density [12]. 
To evaluate the chiral condensate from the eigenvalue density 

Sat (to) = / dxdy — (31) 

1 J x + iy - m 

illustrates the effects of the phase of the fermion determinant in QCD at 
nonzero chemical potential. In particular for p > where the sign 

problem is acute [21], it would be desirable if we could find a real and 
positive reformulation of this integral. 

The Fourier transform of the eigenvalue density along the imaginary 
axis is defined by 

/oo 
dye~ ivt p Nf {x,y), (32) 
-OO 

= /^25(x-x k )e- i «A . 



N 



Since we integrate over all y the Fourier transform is necessarily real. 
The chiral condensate becomes 

p(x,y) 



S(to) = J dxdy 



x + ty — m 

= f dxdy—^ A- / dtj*p{x,t), 

J x + iy - TO 27T J 

Iff f e iyt 

dx / dtp(x, t) I dy- 



2n J J J x + iy — to 

= J ' dx J dtp(x,t)e t{m - x) 9(x-m). (33) 

Now the integrand is real and positive if p is. Below we will see that this is 
the case for the asymptotic limit of the unquenched microscopic density. 
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4.1. The asymptotic limit of the microscopic spectral 
density 

In the limit where the microscopic variables m and /t are much larger than 
unity, the expression for the eigenvalue density simplifies considerably. The 
quenched part of the microscopic eigenvalue distribution is simply constant 
and equal to 1/4/t 2 for Re[i] < 2/t 2 . This leads to a chiral condensate that 
decreases linearly with the quark mass 

= ^ for \m\ < 2/t 2 . (34) 

ZijJi 

In the asymptotic limit the oscillating part of the microscopic eigenvalue 
density reads [12] (here and below to > 0) 

Pu (x,y) = _L_ e (™-*-*y)(™+3*-8A 2 -^)/(8/i 2 )^ ( 35 ) 

Obviously the amplitude grows exponentially with the volume whereas the 
period is inversely proportional to the volume. 

By simple Gaussian integration it follows that the Fourier transform of 
Pu is 

p^t) = ^= e -2A 2 (t-l) 2 -S 2 /2A 2 -™t-(t-2)S+ 2 ™^ (36) 
V27T/2 2 

It is a positive definite and nicely behaved function with maximum at t = 
1 — (to + i)/4/i 2 (see Fig. 3). The contribution to the chiral condensate of 
the unquenched part of the spectral density thus becomes an integral over 
a positive definite function of the real part of the eigenvalue, x, and the 
Fourier coordinate, t, 

H u {m) = ^L= dx dte-^^-^-^/^-^-^eix-m^) 
\/2Txp? J -2ji 2 J -co 

A simple computation yields 

r 2/f 



Eu(m) = — t== / dx6{x - m)e- &2 l 2(i l dqe~ 

\/2TTfl 2 J-2jj? J -co 

= —^[ 2 ' 1 dx8(x-m)e- £2 /^\P^e* 2 / 2 ^ 



2A 2 / 



2[i 2 

dx8(x — to), 

2£ 2 



TO 
2 " 



2/t 



(38) 
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Re[z] 



Fig. 3. Top: The asymptotic limit of the oscillating part of the unquenched eigenvalue 
density. Bottom: The Fourier transform of the asymptotic limit of the oscillating part 
of the unquenched eigenvalue density. In both cases m = 30 and fi = 7. 



The full chiral condensate is obtained by adding the quenched part (34), 
note that the \x dependent term drops out. The message to take away from 
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this derivation is that the calculation of the chiral condensate may simplify 
by using the mixed Fourier transform of the spectral density introduced 
above. 

5. Conclusions 

Chiral symmetry breaking has been linked to complex oscillations on the 
microscopic scale. The cancellation of these violent oscillations, which result 
in the correct value of the chiral condensate, demonstrates the numerical dif- 
ficulties in dealing with the sign problem when the quark mass is inside the 
support of the eigenvalues. At low temperatures this occurs for /i > 
We have shown that the mechanism of these cancellations can be under- 
stood in terms of orthogonality relations for the orthogonal polynomials 
corresponding to a random matrix theory that describes the microscopic 
limit of QCD. The derivation is even valid for finite size random matrices. 
The fact that the oscillations predominantly occur parallel to the imaginary 
axis motivated us to study the Fourier transform of the eigenvalue density 
in the imaginary part of the eigenvalue. In the asymptotic limit the result- 
ing Fourier transformed eigenvalue density is real and positive, and in this 
formulation a probabilistic interpretation of the contributions to the chiral 
condensate is possible. 

It would be most interesting to generalize these results to the p-domain 
of chiral perturbation theory. 
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